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The effect of surfactants on the thermodynamics and kinetics of semiconductor heteroepitaxy 
is briefly discussed. It is argued that the way the surfactants suppress the thermodynamic driving 
force for 3D islanding depends on the mechanism of exchange of overlayer and surfactant atoms. 
If the overlayer atoms occupy bulk-like positions provided by the outgoing surfactant dimers, as is 
the case of Ge/Si(001), large atomic displacements are forcibly inhibited, and the wetting of the 
overlayer by the substrate becomes nearly complete. This complete wetting of the overlayer by the 
substrate means a zero thermodynamic driving force for 3D islanding. Once the thermodynamics do 
not require 3D islanding a change of the growth mode with the temperature should not be observed. 
Thus, the temperature change of the growth mode appears as an indication for the primary role of 
the kinetics, as in the case of Ge/Si(lll). 
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I. INTRODUCTION 



= Vsd + (Eaa - E AB ) = n% D + ^aa * (2) 



The presence of surface active species (surfactants) on 
the surface of growing crystals changes dramatically the 
mode of growth, the rate of nucleation, the transition 
from step flow to two-dimensional (2D) nucleation, the 
structure of steps, the distribution of,stress, surface alloy- 
ing, the concentration of defects, etc.EJ The major effect is 
the suppression of the three-dimensional (3D) islanding 
so that smooth films suitable for fabrication of microelec- 
tronic devices could be produced. Although the initial at- 
tempts to understand the effect of the surfactants on the 
mode of growth were based on thermodynamic grounds,EI 
there is still a debate in the literature between the adher- 
ents of the thermodynamic and the kinetic views. That 
is why we will make an attempt to consider both the 
thermodynamic and the kinetic aspects of the effect of 
surfactants on the thin film growth mode. We will do 
that in the most general way without accounting for the 
particular properties of the materials. We note only that 
in the case of semiconductor growth a complete mono- 
layer of the surfactant is neededE 

II. THERMODYNAMICS 

If we want to change the direction of a process we 
have to change the sign of the thermodynamic force that 
drives it. If we want to prevent a particular process to 
take place we have to suppress the corresponding ther- 
modynamic driving force making it as close to zero as 
possible. The thermodynamic driving force which de- 
termines the occurrence of one or another mechanism of 
epitaxial growth is the difference A/i = fj,(n) — /j!^ d of 
the chemical potential, fi(n), of the overlayer which de- 
pends on the film thickness measured in number n of 
monolayers counted from the interface, the chemical 
potential, fJ% D , of the bulk 3D crystalffl The thickness 
dependence of the film chemical potential fi(n) originates 
mostly from the thickness distribution of the strain due 
to the lattice misfit, but the interaction between the de- 
posit and the substrate, which rapidly decreases with the 
distance from the interface and could be neglected he^ 
yond several monolayers should be also accounted foroEl 
If we deposit a crystal A on the surface of a crystal B in 
absence of a surfactant the thermodynamic driving for.ee 
can be written in terms of surface energies a a and o~gB 

A/i = a 2 [a A + CTAB{n) - <r B ] (1) 

where a 2 is the area occupied by an atom, and o~AB{n) is 
the interfacial energy which includes in itself the misfit 
strain energy and the Attenuation of the energetic influ- 
ence of the substr.ate.El In fact, this is the familiar 3-er 
criterion of Bauer .0 

We can write, the above expression in terms of inter- 
atomic energiesO 



where Eaa and Eab are the energies per atom to disjoin 
a half-crystal A from a like half-crystal A and an unlike 
half-crystal B. It is in fact the adhesion energy Eab 
which includes in itself the thickness distribution of the 
strain energy. In the above equation <Z> = 1 — Eab / Eaa 
is the so called adhesion parameter which accounts for 
the wetting of the overgrowth by the substrate. This is 
the same parameter which enters the work of formation 
of nuclei on a foreign substrate in the classical nucleation 
theory. 

Equations ([!]) and (0) are identical. Eq. can be 
readily obtained from (0) by using the definition of sur- 
face free energy and the relation of Dupre.H The essential 
physics, however, is in Eq. (0J as the surface energies 
are derivatives of the interatomic forces. It shows that 
the thermodynamic driving force for occurrence of one or 
another mode of growth is the relative adhesion/cohesion 
difference, or in other words, the wetting. The incom- 
plete wetting, < $ < 1, is the driving force for 3D 
nucleation and growth, whereas the complete wetting, 
<P < 0, is the driving force for planar growth hy, nucle- 
ation and growth of 2D monolayer height nuclei.El 

When < <P < 1 fi(n) > /J^ D , and 3D islands are 
formed on top of the substrate surface from the very be- 
ginning of deposition. This is the well known Volmcr- 
Weber (VW) growth and fj,(n) goes asymptotically n ta 
/j,® D from above (the curve denoted by VW in Fig. |l|) .□ til 

In the other extreme of complete wetting ( $ < 0) ei- 
ther layer-by-layer (Frank-van der Merwe or FM) growth 
at negligible lattice misfit, or a Stranski-Krastanov (SK) 
growth (layer-by-layer growth followed by 3D islands! at 
a perceptible misfit, are thermodynamically favored. ltEIj 
In the case of FM growth /i(n) goes asymptotically to ^ D 
from below (the curve denoted by FM in Fig. [j]). The 
first monolayer has the lowest chemical potential owing 
to the strongest interaction with the substrate, and a sec- 
ond monolayer can form only after the completion of the 
first one. 

The SK growth mode appears as a result of the inter- 
play of the film - substrate bonding, strain, and surface 
energies. A wetting layer consisting of an integer num- 
ber of monolayers is first formed by a FM mode which 
is driven by a negative thermodynamic driving force A/i 
(complete wetting). Further FM growth becomes unfa- 
vorable as the strain energy accumulates linearly with 
film thickness. In addition, the stronger attraction from 
the substrate, that overcompensates the strain energy, 
disappears beyond several atomic diameters. On top of 
the wetting layer 3D islands form and grow under the 
influence of a positive thermodynamic driving force (in- 
complete wetting). If misfit dislocations (MDs) are in- 
troduced immediately after the completion of the wet- 
ting layer unstrained 3D islands are formed and grow 
on top of it. This is the classical Stranski-Krastanov 
growth in which the work needed to create new surfaces 
is overcompensated by a complete strain relaxation. It 



2 



was found recently that initially coherent (dislocation- 
free) 3D |isJa»ds are formed (coherent Stranski-Krastanov 
growth) ,liirli2l in which the formation of a new surface 
is overcompensated by a gradual strain relaxation. In 
both cases the 3D islands and the wetting layer represent 
necessarily different phases separated by an interphase 
boundary. This boundary is determined by the displace- 
ments of the atoms belonging to the first atomic plane of 
the islands from the bottoms of the corresponding poten- 
tial troughs of the uniformly strained wetting layer. In 
the classical SK swwth the displacements are described 
in terms of MDsJia whereas in the coherent SK growth 
the atoms that are closer to the islands edges are dis- 
placed owing to the islands finite sizeJl3 These displace- 
ments give rise to weaker adhesion, or in other words, to 
incomplete wetting. It is namely this incomplete wetting 
(or weaker adhesion) which appears as the thermody- 
namic driving force for the 3D islanding in the SK mode 
of growth. In other words, we can treat the SK mode as a 
FM mode driven by complete wetting (A/i < 0), followed 
by VW mode driven by incomplete wetting (A/i > 0) (the 
curve denoted by SK in Fig. |l|). 

When considering the effect of the surfactant on the 
mode of growth we note that the chemical potential ^ D 
of the infinitely large crystal does not depend on whether, 
there are impurities adsorbed on its surface or notJiSEZI 
This, however, is not true for sufficiently thin films. If a 
monolayer thick film of A is covered by a monolayer of a 
surfactant (S) atoms the growth mode criterion ([!]) turns 
intoN 

A/i S = a 2 (a AB + a s a - <Jsb) (3) 

where the cr's are the corresponding interfacial energies. 
This expression is equivalent to 

Wr(n) = ML + (E AA - E AB ) - {E SA - E SB ) (4) 

where /^s(n) is the chemical potential of the film covered 
by S atoms. It is immediately seen that if 3D islanding is 
thermodynamically favored in absence of a surfactant, i.e. 
Eaa > Eab, the inequality (^J) may have the opposite 
sign if the adhesion of the surfactant is stronger to the 
overlayer rather than to the substrate and thus the third 
term in ([!]) overcompensates the second one. Then the 
surfactant is " good" from thermodynamic point of view 
if it adheres more strongly to the overlayer, thus changing 
the sign of the thermodynamic driving force. 

This is valid, however, only for the transition from VW 
to FM growth as ^8-^he case of the As mediated growth 
of Si on Ge(001).aO Even in this case it is valid only 
at zero misfit. After several monolayers of A the ener- 
getic contact between the surfactant and the substrate 
B is completely lost (E$b ~> Eg a), an d if the misfit is 
large 3D islands should appear resulting in SK mode. In 
systems that follow the SK growth mode where, beyond 
the wetting layer, A is deposited on strained A^i the 
difference Esa — Es b wiil rbe nearly equal to zgppr-JVIany 
papers both theoreticalc2rE-3 and experimentalE3~E2l have 



been devoted to the problem considering mainly the sur- 
factant growth of Ge on Si(001). However, in most of the 
papers quoted above the mechanism of the exchange of 
the surfactant (As or Sb) atoms with Ge atoms belonging 
to the first monolayer on the Si(001) surface have been 
studied. It is obvious, that such studies cannot explain 
the effect of the surfactant to suppress the 3D island- 
ing as the first monolayer belongs to the wetting layer 
and it will grow in a FM mode pseudomorphous with the 
substrate irrespective of there is a surfactant on top or 
not. 

Obviously, in the SK mode the surfactant should make 
the thermodynamic driving force for 3D islanding on top 
of the wetting layer equal to zero for considerable film 
thickness. As mentioned above in the SK mode the atoms 
which are in contact with the wetting layer should be 
displaced from their " bulk-like" positions either to make 
MDs or near to the islands edges. These displacements 
make the wetting incomplete, or which is the same, make 
the adhesion weaker than the cohesion {Eab < Eaa), 
which is the thermodynamic driving force for 3D island- 
ing. As found in a series of papers the surfactant group-V 
atoms form dimers on Si(Ge)(001) which peside_esactly 
where the Si(Ge) dimers should be locatedJ3tZKin3 Af- 
ter the surfactant and overlayer dimers exchange places 
the latter occupy the epitaxial sites that are provided 
by the outgoing surfactant. The surfactant dimers on 
top do not permit large displacements of the overlayer 
atoms underneath, and in turn incomplete wetting. In 
other words, the surfactant atoms press the overlayer 
atoms into bulk-like positions and do not allow the for- 
mation of misfit dislocations. Thus, the surfactants sup- 
press the thermodynamic driving force for 3D islanding 
by not allowing large atomic displacements. Once the 
thermodynamics do not require 3D islanding a change 
of the growth mode with the temperature should not be 
observed. 



III. KINETICS 

In absence of a surfactant the atoms arriving at the 
crystal surface diffuse on it, join pre-existing steps or is- 
lands, or give rise to new islands. At sufficiently high 
temperatures the diffusivity of the atoms on the ter- 
races is high and they reach the pre-existing steps be- 
fore meeting with each other. This results in a step-flow 
growth. At low temperatures and in turn low diffusivity 
the atoms meet with each other before reaching the steps 
and give rise to 2D nuclei on the terraces. In presence 
of a complete monolayer of S atoms two new phenomena 
take place. First, the surfactant changes the energetics 
of incorporation of atoms into ascending and descending 
steps, and second, the overlayer atoms have to exchange 
places with the S atoms on the terraces or at the steps in 
order to join the crystal lattice. Both phenomena lead to 
a change of the diffusivity of the adatoms across the steps 
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and on terraces, and in turn exert an effect on the direc- 
tion of transport (upward or downward) of atoms. More- 
over, the change of diffusivity on terraces affects strongly 
the kinetics of nucleation and growth of the islands thus 
stimulating either the step flow or 2D nucleation. We 
consider in more detail the effect of the surfactant on the 
kinetics of attachment and detachment of atoms to and 
from the steps, and the exchange-deexchange kinetics on 
terraces. 

Whereas an atom approaching an ascending step joins 
it upon striking, an atom approaching a descending step 
has to lOKeijCome an additional Ehrlich-Schwoebel (ES) 
barrier HJo If the ES barrier is low the overlayer grows 
more or less in a layer-by-layer made. Otherwise, in- 
stabilities of different kind appear S3 Zhang and Lagally 
pointed out the possibility of a reduction of the ES bar- 
rier when overlayer atoms exchange sites with S atoms 
that decorate the steps, in such a way that the last bar- 
rier before descending is shifted dowawards below the 
level of the surface diffusion barrier S3 The latter is in 
good agreement with first principle molecular dynamics 
calculations of displacement of Sb dimer at the step edge 
by a Si dimer through a push-out mechanism.E3 In ad- 
dition, Markov has suggested that atoms that approach 
ascending steps have to overcome an extra energetic bar- 
rier owing to the necessity to displace the S atoms al- 
ready adsorbed at the step edges SB These effects which 
are of purely kinetic origin can suppress the thermody- 
namic driving force for 3D islanding in heteroepitaxial 
growth. 

When thermodynamics predict 3D islanding, (A/i > 0) 
the atoms are more strongly bound to the upper layers 
rather than to the lower layers. In other words thejjrpper 
layers have lower chemical potentials (Fig. |l]).cm This 
results in a gradient of the chemical potential that drives 
the atoms upwards-agd is the driving force for the 2D- 
3D transformation.oEl The appearance of an additional 
barrier at the ascending steps, and the reduction of the 
ES barrier at the descending steps owing to the decora- 
tion of the steps by S atoms, reverse the asymmetry of 
the attachment - detachment kinetics and thus lead to a 
suppression of the thermodynamic driving force for 2D- 
3D transition. A reverse gradient of adatoms of a kinetic 
origin appears that drives the atoms downwards. r^Ucj 
latter results in a planar rather than a 3D growthEj'cHl 
Thus the attachment - detachment kinetics suppress the 
thermodynamic driving force for 3D islanding. At high 
temperature the thermodynamics prevail and SK growth 
takes place irrespective of the presence of the surfactant. 
Surfactant induced FM growth occurs at lower temper- 
atures. Note that this kinetic effect does not require a 
complete monolayer of the surfactant but a small amount 
which is sufficient to decorate the steps. 

In the presence of a complete monolayer of the surfac- 
tant a new kinetic effect appears owing to the necessity 
of the overlayer atoms to join the crystal lattice and of 
the S atoms to float over the surface. These are the phe- 
nomena of exchange and deexchange of overlayer and S 



atoms. Zhang and LagallyEZI first assumed that an ex- 
change process between overlayer and S atoms on the 
terraces should take place on the S precovered surfaces. 
Kandel and Kaxiras assumed later that a deexchange 
process is also possible and can play a significant role 
in epitaxial growth.Ho They calculated the values of 0.8 
and 1.6 eV for the activation energies, E ex and Edex, for 
exchange and deexchange processes, respectively. Ko, 
Chang and Yi computed through first-principles pseu- 
dopotential total-energy calculations the values of the 
barriers for surface diffusion, exchange and deexchange 
for Ge on Ga, As and Sb precovered Si(lll) surface. 
They found that in all cases Edex > E ex but the differ- 
ences are not-as large as the ones calculated by Kandel 
and Kaxiras.E3 

We introduce two time constants which characterize 
the exchange and the deexchange processes both being 
normalized to the time, t\ = 1/F, necessary for depo- 
sition of a completc-JTionolayer, F being the the atom 
arrival frequency.e3'Ej First, this is the mean residence 
time of the atoms on top of the S layer before exchange 
T e x = {F/v)cxp(E ex /kT) where v is the attempt fre- 
quency. Second, this is the mean residence time before 
deexchange of an atom embedded in (or under) the S 
layer r dex = (F/v) exp (E dex /kT). Obviously, if a par- 
ticular time constant is greater than unity the corre- 
sponding process will not occur. As Td ex 3> i~ ex three 
possibilities exist. The first is Tdex *C 1. This extreme 
describes the reversible exchange as the adatoms have 
sufficient time to go back on top of the S layer through 
a deexchange process. A dynamic exchange - deexchange 
equilibrium is established. With the value of Edex =1.6 
eV calculated by Kandel and Kaxiras we find, at 600K, 
Tdex ~ 0.002. Then the exchange - deexchange equilib- 
rium is established in the very beginning of the deposi- 
tion process. e3'E2I The nucleation process takes place in a 
more or less disordered (because of the different size of 
the atoms involved) 2D phase - a monolayer consisting 
of mixed S and overlayer atoms. A considerable frac- 
tion of the incoming atoms remains on top of the S layer 
and diffuse fast to the pre-existing steps. The forma- 
tion, growth and decay of nuclei take place through ex- 
change and deexchange processes. The second extreme 
is T ex « 1 < Tdex- This is the case of irreversible ex- 
change. The incoming atoms rapidly exchange places 
with S atoms and remain buried under the S layer. All 
processes of diffusion, nucleation and incorporation into 
islands and steps occur under (or between) the S atoms. 
The case 1 < r ex <C Tdex should be excluded as it means 
that the S atoms will remain buried under the arriving 
overlayer atoms. 

We calculate further the nucleus density, N$, in both 
cases of reversible and irreversible exchange. This is ex- 
tremely important because of two reasons. First, the 
large number of small islands will promote layer-by-layer 
growth owing to the existeBGe of a critical island size for 
second layer nucleation.BLJ Second, higher nucleation 
rate leads to growth by 2D nuclei, whereas in the oppo- 
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site case step-now growth will take place. In all cases the 
expression for N§ can be written in the form 

N s = N ex V (-?rEs/kT) (5) 

where E$ combines all energy contri butio ns which de- 
pend on the presence of the surfactant ,caE3 and 

N = (^y x exp(*rE a /kT) (6) 

is thje-jpucleus density in absence of a surfactant (Es = 
0)jt3Ei Thus, the island density scales with the ratio v/F, 
the scaling exponent x being a famction of the number 
i of atoms in the critical nucleusJld'Cj In the case of low 
barriers for attachment of atoms to 2D nuclei, the lat- 
ter " feel" (©^ph, other through the diffusion fields around 
them, andNa 

Kandel has recently shown by using a rate equation 
approach that in presence of surfactants which give rise to. 
barriers for incorporation of atoms to the critical nucleic^ 



The scaling exponent (Q) has been independently derived 
by Markov by using a different approach. ej He showed 
that it is characteristic for kinetic regime of growth where 
the adatom concentration between the islands is practi- 
cally constant and the islands " do not feel" the presence 
of each other. ._. 

As has been noted in Ref. (0) the scaling exponent 
(0) varies with i from 1/3 to 1, whereas (||) has values 
larger than unity already at i > 2. Thus, one could dis- 
tinguish between diffusion and kinetic regimes of growth 
if x is smaller or greater than unity. Typical examples 
of the scaling exponent (||) are th&Jiomoepitaxy of Si 
on Sn precovered surface of SLL(111),e3 and of Ge on Pb 
precovered surface of Si(lll)o 

It is immediately seen that the exponential multiplying 
No in Eq. (||) can be smaller or greater than unity de- 
pending on the interplay of the energies involved in Eg. 
If Es < the surfactant will stimulate the nucleus forma- 
tion (Ns ^> Nq). The latter means that a greater density 
of smaller islands will be formed, and the film will grow 
by 2D nucleation rather than by step-flow. What is more 
important is that the greater density of smaller nuclei 
leads to layer-by-layer growth rather than to the mound 
formation. In the opposite case the surfactant will inhibit 
the nucleation and will promote the step-flow growth. 
In the case of reversible exchange Es contains the term 
-i[(E dex -E ex )- AE sd ] where AE sd = E° d - E sd is the 
the difference between the barriers for diffusion on the 
clean substrate, E° d , and on the surfactant passivated 
surface, E sd , (E° d > E sd ). If AE sd overcompensates the 
difference E dex — E ex , we could expect that Es > 0, the 



nucleus formation will be inhibited, (Ns -C No),, and 
surfactant-induced step-flow growth will take place.oEJ 

In the case of irreversible exchange Eq. (||) is still valid 
but the barriers E dex and E ex do not enter the energy 
term Es- The largest term in Es is the increment of 
the surface diffusion barrier —i(E sd — E° d ) where now 
E sd > E° d owing to the fact that the process of diffusion 
under the S layer is inhibited compared with that on the 
clean surface. It is unlikely this term to be overcompen- 
sated by the smaller positive terms, and hence, we can 
always expect that the surfactant will stimulate the 2D 
nucleation rather than the step-flow growth. 

Finally, we discuss briefly the reasons that lead to a 
change of the island density. Usually the latter is at- 
tributed to the change of atoms diffusivity. However, the 
change the island density might be due to other reasons. 
The nucleation rate depends on the Gibbs free energy for 
nucleus formation which is in turn a function of the nu- 
cleus size. The presence of the surfactant could lead tci-a 
change of the number of atoms in the critical nucleus.Eil 
As a result the values of Es and E in Eqs. (||) and (||), 
as well as the scaling exponent x wm change in one and 
the same temperature interval. In addition, the regime of 
growth can change from diffusion to kinetic, so that the 
scaling exponent (fy will be replaced by (||) . In this case, 
the island density will increase by orders of magnitude 
even if i remains constant. 



IV. GROWTH OF GE ON SI 

We discuss in more detail a "model!, system such 
as Ge/Si. As pointed out by KaxirasJ^Zl one and the 
same surfactant can have qualitatively different effect de- 
pending on the crystallographic orientation of the sub- 
strate. On Si or Ge(001) the surfactants forro_di«icrs 
which do not break during the exchange process .EirEZJ On 
(111) surface three different geometries (substitutional, 
trimers, and zig-zag chains) of the group-V adsorbates 
(P, As, Sb) are possible.Ell In the first lxl structure one 
monolayer of S atoms is incorporated into the upper part 
of the (111) bilayer, whereas in the last two structures a 
monolayer of S atoms is bonded on top of a complete 
(111) bilayer. In addition, the growth of the two sur- 
faces of Si or Ge differs drastically—Whereas the (001) 
surface grows by single monolayers, E3 the (111) surfaces 
grow by bilayers.Ea Thus, we have completely different 
mechanisms of exchange of S and overlayer atoms. In 
the case of (001) surfaces the growth unit is a dimer so 
that an exchange of S and ovepiayer dimers takes place 
through a single monolayer .EiTED On (111) surfaces the 
mechanism of exchange depends on the surface geometry 
of the surfactant. For example, if the surfactant geometry 
is substitutional we could suppose that single overlayer 
atoms should displace single S atoms from their positions 
in the upper monolayer that constitutes the bilayer. Af- 
ter the complete building of the upper monolayer atoms 
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belonging to the next lower monolayer should be incor- 
porated first in order to serve as a template for the S 
atoms in the next upper monolayer. Thus one surfactant 
monolayer should be displaced by a complete bilayer of 
the overgrowth the process taking place by exchange of 
single atoms. Comparing both surfaces we could specu- 
late that on (111) surface the surfactant could not sup- 
press successfully the thermodynamic driving force for 
3D islanding by not allowing large atomic displacements 
(formation of misfit dislocations). 

We consider first the growth of Ge on Si(lll). Hwang, 
Chang and Tsong studied this system by using Pb as a 
surfactant. They established that the Pb atoms make 
a substitutional lxl structure on top of the Ge sur- 
face. It was found that the saturation nucleus density in 
the submonolaiier regime scales with the deposition flux 
with x = 1.76x3 This clearly shows that in the Pb me- 
diated epitaxy of Ge on Si(lll) the scaling exponent is 
given by Eq. (|^) and the growth proceeds in a kinetic 
regime. In other words, significant energy barriers ap- 
pear for attachment of atoms to the critical nuclei owing 
to the necessity of exchange of overlayer and S atoms. 
The same scaling exponent (x = 1.8) has been obtained 
in the homoepitaxial growth of Si on Sn precovered sur- 
face of Si(lll)^3 after comparing the experimental data 
for the critical terrace wid^h-jfor step flow growth with 
the corresponding theoryM^'Ed Under clean conditions in 
the same system Voigtlander et al. obtained a scaling 
exponent smaller than unity (x = 0.85) J£j 

The observations of Hwang et al. show unambiguously 
that Pb has a prominent kinetin. effect on the submono- 
layer growth of Ge on Si(lll)o In a further study the 
same authors found that 3D islanding is suppressed be- 
low 470K. At temperatures higher than 470K 3D island- 
ing takes place.Eil This observation could be explained by 
the suppression of the thermodynamic driving force for 
3D islanding by the attachment - detachment kinetics 
mentioned above. The upward surface flux is hindered 
whereas the downward flux is enhanced. We could con- 
clude that the suppression of 3D islanding in Ge growth 
on Pb precovered Si(lll) is due more to kinetic rather 
than to thermodynamic reasons. 

Voigtlander and Zinner-studied the growth of Ge on 
Si(lll) mediated by Sb.E3 They established the same 
substitutional structure of the Sb atoms on the surface 
of Ge. They found that a transition from surfactant me- 
diated layer-by-layer growth to the equilibrium SK mode 
takes place above twice higher temperatures (900K) as 
compared with the Pb mediated growth. The higher tem- 
peratures could be related to stronger chemical bonding 
but qualitatively the surfactant effect is the same. _ 

The system Ge/Si(001) is the most studied one.0 Ide 
has observed that after annealing at about 800K for 10 
min the 2D Ge islands formed in submonolayer regime on 
As precovered Si(001) did not change whereaslhe islands 
deposited under clean conditions disappeared^ This ob- 
servation shows that 2D islands formed under the surfac- 
tant layer are more stable than those without S atoms 



on top of them. Antimony mediated deposition and an- 
nealing experiments carried out by Osten et al. showed 
that Sb prevents 3D isJanding, but alsq_can smooth an 
already islanded filmB3 Katayama et alx3 found that as 
the coverage of the surfactant increases, the intermixing 
of Ge and Si, as well as the nucleation and the growth 
of macroscopic Ge islands, are suppressed. Both obser- 
vations could be explained with the forced occupation of 
epitaxial sites by the Ge atoms under the influence of the 
S dimers on top of them. On the other hand, Ide also es- 
tablished that at 800K the Ge atoms reach and join the 
preexisting steps of Si (step flow growth) on the clean 
surface- whereas 2D islands are formed in the presence 
of As.Ea The latter unambiguously demonstrates the ki- 
netic effect of the surfactant on the mode of growth by 
decreasing the diffusivity of the incoming atoms. 

Most of the studies of surfactant mediated growth 
of Ge or .Si fnoi ) siirfapp are carried out at or near to 
500°CByH!3'HHl£3 To the author's knowledge higher 
temperature deposition has been carried out by Tromn 
and Reuter (630°C),B and by Horn von Hoegen et al.H 
at 700°C. No 3D islanding except for mound formation 
has been noticed on top, of the wetting layer at the high- 
est temperature used.EJ A change of the growth mode 
with temperature as in the case of Ge/Si(lll) has not 
been reported so far£3 This can be interpreted as fol- 
lows. The surfactant suppresses successfully the thermo- 
dynamic driving force for 3D islanding by not allowing 
large atomic displacements and the growth continues in 
a layer-by-layer mode (with the unavoidable roughness 
to relieve the strain) untill MDs are introduced at the 
interface. Once the thermodynamics do not require 3D 
islanding a change of the growth mode with the temper- 
ature should not be observed. 

Jenkins and Srivastava carried out first principle den- 
sity functional theory (DFT) calculations-^ the mode 
of growth of Ge on Sb precovered Si(001).E-3 They stud- 
ied the behavior of the chemical potential of the Ge film 
relative to the bulk chemical potential, /J^ Dl in absence 
and presence of Sb as a surfactant. They found that in 
absence of Sb the chemical potentials of the first three 
monolayers of Ge, which belong to the wetting layer, are 
smaller (more negative) than /i° D (complete wetting). 
The chemical potential of the forth monolayer becomes 
more positive than fi® D (incomplete wetting) which is an 
indication of 3D islanding after the first three monolayers. 
In presence of Sb the chemical potential of the first mono- 
layer is slightly more negative that [i® D but the chemical 
potentials of the remaining three monolayers are equal 
to /ig D within the accuracy of the calculations. An er- 
ror is introduced in the calculations by not allowing the 
Ge overlayers to relax laterally. Thus after the exhange 
of Ge and Sb dimers the Ge atoms occupy the bulk-like 
positions left by the outgoing Sb dimers and the surfac- 
tant successfully suppressed the thermodynamic driving 
force for 3D islanding. In fact, the surfactant replaced 
the curve SK by the curve FM in Fig. (|I|). 

Similar first principle DFT calculations have been car- 
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ried out by Gonza- les-Mendez and Takeuchi for the. 
growth of Si on Ge(OOl) by using As as a surfactant.Ej 
In absence of a surfactant Si grows-Ln VW mode as 3D 
islands directly on the Ge surfaceH In presence of As 
the overlayer prefers to grow in a layer-by-layer mode 
due to thermodynamic reasons, the energy of the 3D is- 
lands being always lower than that of the complete mono- 
layers. The deposit atoms are almost in bulk-like posi- 
tions. Equally important, the authors established that 
As greatly reduces the Si and Ge intermixing. This is 
another indication of the thermodynamic influence of the 
surfactant as the alloying also requires large atomic dis- 
placements. 



the mechanism of exchange of overlayer and surfactant 
atoms. Undoubtedly the third element changes the sur- 
face energies but it is the change of the relative adhesion 
of the overlayer to the substrate that is crucial from ther- 
modynamic point of view. 
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V. CONCLUSION 



The presence of a third element in the system unavoid- 
ably changes the kinetics of growth. Hence, we have 
to decide whether the kinetics is solely responsible for 
suppression of the 3D islanding, or the thermodynam- 
ics plays a decisive role, with the kinetics exerting an 
additional effect. We argue that this question could be 
answered by studying the temperature dependence of the 
capability of the surfactant to suppress 3D islanding. If 
the kinetics play a decisive role we should expect a change 
of the mode of growth from layer-by-layer at lower tem- 
perature to 3D (SK or VW) islands at higher temper- 
ature. The physical reason is that at high temperature 
the system is closer to equilibrium whereas at lower tem- 
perature the kinetic effects suppress the thermodynamic 
driving force for 3D islanding. This is most probably 
the case of the system Ge/Si(lll) where one monolayer 
of the surfactant should be displaced by two monolayers 
of the overgrowth and the exchange process takes ran- 
domly place by single atoms. The S atoms do not press 
the overlayer atoms into bulk-like positions and forma- 
tion of misfit dislocations or displacements of the edge 
atoms is not inhibited. We have seen that the growth 
is accompanied by typical kinetic phenomena as the ki- 
netic regime of growth (x > !)■ O n the contrary, the 
system Ge/Si(001) demonstrates the decisive role of the 
thermodynamics. The kinetics play an auxiliary role in 
the same direction. The physical reason for suppression 
of the 3D islanding is that a monolayer of the surfactant 
is displaced by a monolayer of the film during growth, 
and the exchange process occurs by dimers rather than 
by single atoms. As a result the overlayer atoms occupy 
nearly bulk-like positions provided by the outgoing sur- 
factant dimers. Large atomic displacements are forcibly 
prohibited which in turn leads to complete wetting and a 
zero thermodynamic driving force for 3D islanding. It is 
thus clear that the third element in the system does not 
act as a surfactant at all. 

We conclude that the relative weight of the effect of 
the third element on the thermodynamics and kinet- 
ics in semiconductor heteroepitaxy depends mostly on 
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FIG. 1. Schematic dependence of the film chemical po- 
tent- ial on the film thickness in number of monolayers for 
the three modes of growth denoted at each curve: VW - 
Volmer- Weber, SK - Stranski-Krastanov, and FM - Frank-van 
der Merwe. 
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